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Major Drivers of Current Work 
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Impact on component design => Rapid turn-around of high-fidelity analysis 
Increase durability /safety => Accurate quantification of flow 

(i.e. prediction of flow-induced vibration) 

Impact on system performance => More complete systems analysis 

using high-fidelity tools 
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Current Challenges 
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- Time-integration scheme, convergence 

- Moving grid system, zonal connectivity 

- Parallel coding and scalability 
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- Scalability (algorithmic & implementation) 

- Portability, transparent coding, etc. 

Co m p u t i n q res o u r c 6 s 

- "Grid" computing will provide new computing resources for 
problem solving environment 

- High-fidelity flow analysis is likely to be performed using "super 
node" which is largely based on parallel architecture 
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INSSb - Incompressible A/-5 Solver 





** Parallel version ; 

•MPI and MLP parallel versions 

• Structured, overset grid orientation 

• Moving grid capability 

• Based on method of artificial compressibility 

• Both steady-state and time-accurate formulations 

• 3 rd and 5 th -order flux difference splitting for convective terms 

• Central differencing for viscous terms 

• One- and two-equations turbulence models 

• Several linear solvers : GMRES, GS line-relaxation, LU-SGS, 

G5 point relaxation, ILU(O) 


1982-1987 Original version of INS3D - Kwak, Chang 
1988-1999 Three different versions were devoped : 

INS3D-UP / Rogers, Kiris, Kwak 
INS3D-LU / Yoon, Kwak 
INS3D-FS / Rosenfeld, Kiris, Kwak 
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Time Accurate Formulation 



• Time-integration scheme 









• Introduce a pseudo-time level and artif icial compressibility 
• Iterate the equations in pseudo-time for each time step until 
incompressibility condition is satisfied. 


• Solve auxiliary velocity field first, then enforce 
incompressibility condition by solving a Poisson equation 
for pressure. 
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Space Shuttle Main Engine Turbopump 
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Inlet Guide Vane (IGV) 

•15 Blades 

•Pitch, p = 24 degrees 

•Blade Inlet Angle (mean), p iGV 4 = 90 degrees 
•Blade Exit Angle (mean), P IGV 2 = 45 degrees 

Clearance between IGV and Impeller, x = 0.12 inches 

Impeller 

•6+6+12 Unshrouded Design 
•Pitch, p = 60 degrees 

•Blade Inlet Angle (mean), p impl = 23 degrees 
•Blade Exit Angle (mean), p imp 2 = 65 degrees 
•Clearance between blade LE and Shroud, r = 0.0056 inches 
•Clearance between blade TE and Shroud, x = 0.0912 inches 

Clearance between Impeller and Diffuser, r = 0.050 inches 

Diffuser 
•23 Blades 

•Pitch, p = 15.652 degrees 

•Blade Inlet Angle (mean), p difl = 12 degrees 

•Blade Exit Angle (mean), p dif 2 = 43 degrees 
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INS30 Parallelization 


19.2M Points 

MPI-OpenMP Hybrid 
A- NAS-MLP 


19.2M Points 

MPI-OpenMP Hybrid 
A- NAS-MLP 


12 Groups 


20 Groups 


Number of CPUs 


Number of CPUs 
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Inlet Guide Vane 
15 Blades 
23 Zones 
6.5 M Points 


Diffuser 
23 Blades 
31 Zones 
8.6 M Poi 



60 Zones / 19.2 Million Grid Points 
Overset connectivity : DCF (B. Meakin) 
Less #»an 156 orphan points. 






Scripting Capability 





SCRIPTING CAPABILITY FOR GRID GENERATION 


> Require expertise to build scripts the first time 

> Allow rapid re-run of entire grid generation process 

> Easy to do grid refinement and parameter studies 

> Easy to try different gridding strategies 

> Documentation of gridding procedure 

> Written in Tel scripting language 

> works on UNIX, LINUX and WINDOWS 

> integer and floating point arithmetic capability 

> modular procedure calls 

> easy to add GUI later if needed 


17 


Scripting Capability 


INPUT AND OUTPUT 

Current example: one script for each component 

(IGV, Impeller and Diffuser) 



> profile curve for hub and shroud in PLOT3D format 
(rotated by script to form surface of revolution) 

> blade and tip surfaces in PLOT3D format 

> Parameters that can be changed 

- global surface grid spacing (on smooth part of geometry) 

- local surface grid spacing (leading/trailing edges, etc.) 

- normal wall grid spacing (viscous, wall function) 

- marching distance 

- grid stretching ratio 

- number of blades 




utput 

> overset surface and volume grids for hub, shroud, blades 
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Scripting Capability 




INLET GUIDE VANES AND DIFFUSER 


Old IGV New IGV Old DIFF New DIFF 


No. of points (million) 7.1 


Time to build 


1/2 day 10 sec. 1/2 day 8 sec. 


Script timings on new grids based on SGI R12k 300MHz processor 
Time to build script = 1 day for IGV, 1 day for DIFF 
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Scripting Capability 



FUTURE PLANS FOR SCRIPTING 


> Complete domain connectivity capability in scripts 
(X-ray maps and DCF input file creation) 

> Flow solver input creation in scripts 

> Perform more tests on different parameters 


> Perform tests on different geometries, e.g., volute, inducer 

> Improve robustness (error traps, wider range of cases) 


> Generic template for each component 

> Graphical interface front end 
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VELOCITY MAGNITUDE 




rotations are completed. 

*One physical time-step requires 
less then 20 minutes wall time 
with 128 CPU's on SGI Origin 
platforms. One complete rotation 
requires one-week wall time. 

Code optimization is currently 
underway. For small case, 50Yo 
improvement is obtained by 
employing a better cash usage in 
the code. Less than 10 minutes 
per time step will be obtained by 
the end of September 2001. 
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STATIC/DYNAMIC STRESS ANALYSIS 
FOR TURBOPUMP SUB-SYSTEMS 
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FLUID/STRUCTURE INTERFACE 



• LUMPED LOAD APPROACH 

- FAST, NEEDS FINE GRIDS, ADEQUATE FOR UNCOUPLED METHOD 


• CONSISTENT LOAD APPROACH (CONSERVES LOADS) 
- ACCURATE FOR COUPLED METHODS, EXPENSIVE 




CONSISTENT LOAD APPROACH USING VIRTUAL SURFACE VALIDATED IN ENSAERO 



Gy 


ruswamy 
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STRUCTURES 






STRUCTURES WILL BE MODELED USING 

BEAM, PLATE, SHELL AND SOLID FINITE ELEMENTS 


COARSE GRID 
230 NODES 
414 FE 

1196 DOF k 


TYPICAL 
STRUCTURAL 
MODE AT 12KHz 
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Summary 



•Unsteady flow simulations for RLV 2 nd Gen baseline turbopump for one and 
half impeller rotations are completed by using 34.3 Million grid points model. 

•MLP shared memory parallelism has been implemented in INS3D, and 
benchmarked. Code optimization for cash based platforms will be completed 
by the end of September 2001. 

•Moving boundary capability is obtained by using DCF module. 

• Scripting capability from CAD geometry to solution is developed. 

• Data compression is applied to reduce data size in post processing. 

• Fluid/Structure coupling is initiated. 
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